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Abstract: Topography plays an important role in determining the glacier changes. However, topography 
has often been oversimplified in the studies of the glacier changes. No systematic studies have been 
conducted to evaluate the relationship between the glacier changes and topographic features. The present 
study provided a detailed insight into the changes in the two branches (east branch and west branch) of 
Urumqi Glacier No. 1 in the Chinese Tianshan Mountains since 1993 and systematically discussed the 
effect of topography on the glacier parameters. This study analyzed comprehensive recently observed data 
(from 1992/1993 to 2018/2019), including mass balance, ice thickness, surface elevation, ice velocity, 
terminus, and area, and then determined the differences in the changes of the two branches and explored 
the effect of topography on the glacier changes. We also applied a topographic solar radiation model to 
analyze the influence of topography on the incoming shortwave radiation (SWin) across the entire glacier, 
focusing on the difference in the SWin between the two branches. The glacier mass balance of the east 
branch was more negative than that of the west branch from 1992/1993 to 2018/2019, and this was 
mainly attributed to the lower average altitude of the east branch. Compared with the west branch, the 
decrease rate of the ice velocity was lower in the east branch owing to its relatively increased slope. The 


narrow shape of the west branch and its southeast aspect in the earlier period resulted in a larger glacier 
terminus retreat of the west branch. The spatial variability of the SWin across the glacier surface became 
much larger as altitude increased. The SWin received by the east branch was slightly larger than that 
received by the west branch, and the northern aspect could receive more SWin, leading to glacier melting. 
In the future, the difference of the glacier changes between the two branches will continue to exist due to 
their topographic differences. This work is fundamental to understanding how topographic features affect 
the glacier changes, and provides information for building different types of relationship between the 
glacier area and ice volume to promote further studies on the basin-scale glacier classification. 
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1 Introduction 


As large solid reservoirs, mountain glaciers have a vital impact on regional water resources, 
ecosystems, and economic societies (Zemp et al., 2015; RGI Consortium, 2017; IPCC, 2019). 
Over the last 50 a, climate warming has caused most glaciers to shrink on a global scale, and this 
shrinkage has been accelerated since the end of the last century (Vaughan et al., 2013; Farinotti et 
al., 2015; Zemp et al., 2015). The glacier mass loss in the Tianshan Mountains is mostly 
unsustainable, and ice volume in this region will reduce by at least 30.00% by 2100 (Miles et al., 
2021). Nevertheless, glacier changes are not homogeneous worldwide (Yao et al., 2012; Vaughan 
et al., 2013; Zemp et al., 2015, 2019), and even in the same region under the same regional 
climate conditions, the changes in glaciers may differ because of differences in their dimensions, 
shape types, and topographic characteristics (Wang et al., 2020). For example, glaciers in the 
Tomor Peak of the Tianshan Mountains have been investigated in depth; specifically, glaciers 
have demonstrated quite different shrinkage rates, with one terminus even becoming almost stable 
owing to differences in topography and debris coverage. Therefore, it is important to understand 
the effect of topography on the changes in glaciers. 

The surface topographic settings in which most small glaciers (<0.400 km?) are located could 
result in a very limited glacier area, suggesting that topographic settings lead to the protection of 
glaciers by direct solar radiation and the settings are favourable for glacier preservation (DeBeer 
and Sharp, 2007). Additionally, some glaciers with relatively high altitude caused by topographic 
settings do not extend to low altitude, which reduces glacier mass loss and benefits glacier 
preservation (DeBeer and Sharp, 2007, 2009). Thus, the local surface topographic settings of 
small glaciers have a significant effect on their responses to climate change. The basal topography 
also influences the retreat rate of glaciers. For example, basal topography strongly modulates the 
Humboldt Glacier in northern Greenland, which can produce substantial variations in the glacier 
retreat rate at seasonal to decadal timescales (Carr et al., 2015). Topography significantly alters 
the incoming shortwave radiation (SWin) received at the surface of valley glaciers in High 
Mountain Asia; during the summer months (May—August), the net shortwave radiation of a 
glacier is one of the main components of the surface energy balance, usually accounting for 
84.00% or more of available energy on the surface (Wang et al., 2020). Thus, changes in the 
amount of the SWin on the surface could alter the overall radiation and consequently have a 
significant influence on the surface energy balance. The intensity of the SWin received on the 
surface of a glacier is primarily a function of latitude and the time of year, with components such 
as topographic shading, slope, and aspect controlling the distribution of radiation at a local scale. 
Therefore, it is an important factor not only contributing to the surface energy balance, but also 
influencing glacier response and mass balance estimates (Olson and Rupper, 2019). 

In the middle part of the Chinese Tianshan Mountains, Glacier No. 1 at the head of the Urumqi 
River (namely Urumqi Glacier No. 1) has been monitored since 1959, and it was found that the 
glacier was completely separated into two branches in 1993. The observation results showed that 
the two branches also have different behaviours in regards to how they change. Although many 
studies have reported observations and study results for this glacier (Han et al., 2006; Jing et al., 
2006; Li et al., 2010; Xu et al., 2019), the effect of topography on the changes of this glacier has 
not been addressed in depth (Wang et al., 2018). Therefore, in this study, we aim to present 
comprehensive recently observed data (including mass balance, ice thickness, surface elevation, 
ice velocity, terminus, and area) to discuss the effect of topography on the changes of Urumqi 
Glacier No. 1 based on the differences of the two branches. Importantly, the spatial heterogeneity 
of the SWin across topography, aspect, and slope was also assessed to accurately discuss the 
distribution of point-source SWin observations across the entire glacier. This study evaluated the 
impact of topographic features on the glacier changes in combination with climate change, so as 
to make the cognitive process of glacier changes more comprehensive. 
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2 Study area 


Urumqi Glacier No. 1 (43°07'N, 86°49’E), located in the central region of the Chinese Tianshan 
Mountains in Central Asia, is a cirque-valley glacier with two branches (Fig. 1). This glacier was 
first investigated in 1959 because it was continuously monitored by the Tianshan Glaciological 
Station. Its area and length were 1.910 km? and 2.20 km, respectively, in the early 1960s, based on 
the first mapping in 1962. Since then, it has experienced continuous shrinkage and was completely 
separated into two independent branches in 1993, presently covering an area of 1.521 km’, of 
which the east and west branches occupied 0.970 and 0.551 km’, respectively, based on 
observations in 2018. 


West branch East branch 
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Fig.1 Location of Urumqi Glacier No. 1 (a), overview of the two branches of Urumqi Glacier No. 1 (b), and 
layout of the surface topographic characteristics of the west branch (c) and east branch (d) 


The glacier region is mainly controlled by three dynamic factors: a westerly jet in the upper 
troposphere, the Siberian anticyclonic circulation, and the westerly cyclonic disturbance. 
According to the observation data of the Daxigou Meteorological Station (3539 m a.s.1.), which is 
3 km southeast of Urumqi Glacier No. 1, the mean annual temperature was approximately —4.9°C 
and the annual precipitation was 473.0 mm during the period of 1959-2019. The mean annual 
temperature increased at a rate of 0.24°C/10a, and the annual precipitation increased gradually at 
a rate of 21.8 mm/10a, based on the linear analyses shown in Figure 2. The main ablation period 
is from May to August, and precipitation mainly originates from the water vapour carried by the 
westerly wind. The precipitation from May to August accounts for 78.00% of the annual total 
amount, with the main types of solid fall, such as snow, hail, and sleet. 


3 Observation data and processing method 


The two branches of Urumqi Glacier No. 1 were observed in terms of their mass balance, ice 
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Fig.2 Changes in the annual average temperature and annual precipitation at the Daxigou Meteorological 
Station from 1959 to 2019 


the observation items mentioned and the data of those items. Table 1 provides all the observation 
items mentioned and the data of those items. As Urumqi Glacier No. 1 was separated in 1993, we 
extended some data items, such as the terminus, before 1993 to compare the situations before and 
after the separation of the glacier. In this paper, we briefly introduced the observations, data 
processing, and related uncertainties in recent years. 


Table 1 Data of Urumqi Glacier No. 1 used in this study 


Item Method/Source Measurement frequency Unit Data period 

Mass balance Stake/snow pit Once a month in ablation season m w.e. 1993-2019 

Ice thickness Ground-penetrating radar Once every 5—6 a m 2001, 2006, and 2012 
Stereo-photography, theodolite, 
total station, and RTK GPS Once every few years - 1994-2012 

DEM Terrestrial laser scanning Every year 2015-2017 
Unmanned aerial vehicle End of August 2018 ia 2018 
Theodolite 1994-2001 

Ice velocity Total station Every year m/a 2002-2006 
RTK-GPS 2007-2019 

Terminus Tape and RTK-GPS BE ani eado the ablation m/a 1980-2020 
Stereo-photography, theodolite, 
total station, and RTK GPS Once every few years 1994-2012 

E Terrestrial laser scanning Every year km? 2015-2017 
Unmanned aerial vehicle End of August 2018 2018 

SWin CNR4 radiation sensors Every hour W/m? 30 April r April 


Note: m w.e., meter water equivalent; DEM, digital elevation model; RTK, Real-Time Kinematic; GPS, Global Positioning System; 


SWin, incoming shortwave radiation; -, dimensionless. 
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3.1 Observation data 
3.1.1 Mass balance 


The mass balance of Urumqi Glacier No. 1 has been conventionally observed using the 
stake/snow pit method during the ablation season from 1 May to 31 August since 1959. The mass 
balance from 1993 to 2019 was used for the subsequent analyses in this study. The number of 
stakes differed in various years, with a total of 42 or 43 stakes distributed across the entire glacier 
in recent years (Fig. 1b). The observation items for each mass balance stake included the vertical 
height of the stakes over the glacier surface and the stratigraphy of the snow pits, including the 
thickness and density of the firn layer and the thickness of the superimposed ice. The density of 
the superimposed ice was assumed to be 900 kg/m? owing to the difficulty of the in situ 
measurements (Xu et al., 2019). 

3.1.2 Ice thickness 

Integrated observations using the ground penetrating radar (GPR) equipment were carried out to 
determine the ice thickness of Urumqi Glacier No. 1 in 2001, 2006, and 2012. The GPR 
equipment comprised an enhancement radar system (pulse EKKOPRO 100 A, Sensors and 
Software Inc., Mississauga, Canada) equipped with 100 MHz resistively loaded dipole antennae. 
The field survey was conducted with the antennae oriented perpendicular to the survey lines, 
ensuring that the maximum power was transmitted in-line. The reflections from the offline were 
minimal, although it was possible for high-conductivity surface features to produce reflections 
(e.g., survey poles and data loggers) (Sensors and Software Inc., 1994; Murray et al., 1997). The 
positions of the GPR survey points were concurrently determined using the Real-Time Kinematic 
(RTK) Global Positioning System (GPS) (Unistrong E650, Beijing UniStrong Science and 
Technology Co. Ltd., Beijing, China). The GPR measurement paths included seven transverse 
profiles along the stake rows and one longitudinal profile along the centre lines in the east and 
west branches, respectively. The GPR data were processed using the EKKO View Deluxe 
software (professional version developed by the Sensors and Software Inc., Mississauga, 
Canada), and the corresponding two-dimensional radar images of the survey paths were obtained. 


3.1.3 Surface elevation 


The surface elevation was obtained using a theodolite in 2001, a total station in 2006, and an RTK 
GPS in 2009 and 2012. Specifically, the Unistrong E650 GPS (produced by Beijing UniStrong 
Science and Technology Co. Ltd., Beijing, China) was used in 2009 and 2012. Digital elevation 
models (DEMs) with a resolution of 1 m were established by interpolating GPS survey data. A 
map from the terrestrial stereo-photographic survey of 1994 was also used to derive the DEM. 
From 2015 to 2017, a terrestrial laser scanning system, Riegal VZ®-6000 (Riegal Company, 
Horn, Austria), was used to survey the glacier surface elevation from four scan positions at the 
glacier terminus every year (Xu et al., 2019). An unmanned aerial vehicle (UAV; MATRICE 200, 
DJ-Innovations, Shenzhen, China) survey was conducted at the end of August 2018, covering 
nearly the entire Urumqi Glacier No. 1. Then, we used the Pix4Dmapper software (Pix4D S.A., 
Prilly, Switzerland) to generate a glacier DEM in 2018 with a spatial resolution of 1 m, according 
to the method of Li et al. (2022). The DEMs were resampled at the same spatial resolution. We 
further compared the DEMs of glaciers at different periods to determine the glacier surface 
elevation change over the respective time periods. According to Nuth and Kaab (2011), it is 
necessary to co-register the DEMs before extracting the glacier surface elevation changes. 

3.1.4 Ice velocity 

Traditionally, the ice velocity can be obtained by installing a network of stakes and surveying the 
changes in the positions of the stakes on the glacier surface. Ablation stakes were also used for the 
ice velocity measurements for this glacier. Previously, theodolite and total station instruments 
were used, whereas GPS did not come into use until 2007. The GPS was then used to measure the 
stake positions after 2007. The measurement methods varied based on the year, as explained by 
Wang et al. (2018). The GPS survey method was described in Section 3.1.3. We calculated the 
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displacement vectors in the X and Y directions of the ablation stakes based on two measurements 
taken within a certain period. The stake displacement was then considered as the ice velocity at 
the corresponding stake position. Similar to the other glacier observed items mentioned above, the 
ice velocity change during the period of 1994—2019 was analyzed. 


3.1.5 Terminus position and glacier area 


Changes in the terminal position of Urumqi Glacier No. 1 were determined using both the GPS 
survey and tape measurements. Several sites around the terminus edge were measured using the 
GPS at the beginning (1 May) and end (31 August) of the ablation season, and the distance of 
each site to a stationary landmark was simultaneously measured using a tape measure. These 
actions were taken for both the east and west branches. 

The glacier area was determined by terrestrial stereo-photographic mapping in 1994 and 2000, 
theodolite in 2001, total station in 2006, and RTK GPS in 2009 and 2012. The measurement 
procedures using theodolite and total station instruments have been introduced in previous studies 
(Dahl-Jensen et al., 1986; Wang et al., 2014). Since 2015, laser-scanning systems have been used 
to determine the glacier area (Xu et al., 2019). In 2018, the glacier area was surveyed using a 
UAV. 

3.1.6 Topographic characteristics 

To analyze the effect of topography on the changes of glacier, we derived topographic 
characteristics including surface elevation distribution, slope, and aspect from the DEMs in 1994 
and 2018 and processed them by the 3D Analyst Tools in the ArcMap 10.3 software (Environment 
System Research Institute, ESRI, Redlands, USA). The subglacial topography was constructed by 
subtracting the ice thickness surveyed in 2012 from the glacier surface elevation during the same 
period. The basal topography was analyzed by assuming that the subglacial topography remained 
unchanged during these periods. 

3.1.7 Incoming shortwave radiation (SWin) 

To investigate the distribution of the SWin across the entire glacier, we observed a point-source 


SWin at an automatic weather station (4025 m a.s.1.) deployed on Urumqi Glacier No. 1 from 30 
April 2018 to 30 April 2019. 


3.2 Processing methods 


3.2.1 Mass balance change 


The mass balance of a single point (b) can be calculated using the method of Paterson (1994) as 
follows: 


b=b, +b,; +b,; , (1) 


where bs, bei, and bs; are the mass balances of the snow, glacier ice, and superimposed ice, 
respectively, with the unit of m w.e. (meter water equivalent). 

Based on point observations, we can calculate the specific mass balance using one of the 
following two methods. One method is to interpolate the point mass balance by drawing equal 
mass-balance contour lines. The specific mass balance (bn) is obtained as follows: 


b, = IS sb, (2) 
S 


where s is the total area of the glacier (km7); s; refers to the area between the two contour lines 
(km7); and b; is the corresponding mass balance (m w.e.). 

The second method is based on the repeated point measurements at different altitudes. The 
specific mass balance (b,) can be calculated by extrapolating the point mass balance to the entire 
glacier as a linear function of altitude: 


b= > sib! /s, (3) 
i=l 


where n is the total number of altitude bands; s’ 


i 


is the area of the altitude band (km?); and b; 
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is the corresponding mass balance (m w.e.). 

3.2.2 Ice thickness change 

The ice thickness (h; m) can be calculated using the two-way travel time (ts, with the unit of ns 
(nanosecond)) of the radar wave from the vertical axis and the velocity of the radar signal in 
glacier (v; m/ns), as shown in Equation 4: 


t 
h==xv. 4 
> (4) 


The velocity of the radar signal in mountain glaciers was in the range of 0.167—0.171 m/ns 
(Glen and Paren, 1975; Robin, 1975; Narod and Clarke, 1994). In this study, we set the velocity 
value of 0.169 m/ns based on field tests. 

3.2.3 Surface elevation changes 

To determine the glacier surface elevation changes over the respective time periods, we compared 
the DEMs of Urumqi Glacier No. 1 at different periods, as mentioned in Section 3.1.3. It is 
necessary to co-register the DEMs before extracting the glacier surface elevation changes (Nuth 
and Kaab, 2011). The uncertainties of the DEMs are usually estimated in two ways: (1) 
comparing the DEM values with an independent set of GPS points in the non-glaciated area; and 
(2) simultaneously rectifying the horizontal displacements and vertical biases based on the 
cosinusoidal relationship between the surface elevation difference and topographical parameters 
(slope and aspect). 

To analyze the change in the ice thickness during the period of 1994-2018, we attempted to 
synthesize data from different sources, including changes in thickness and surface elevation. 
Finally, the relationship between the ice thickness change and surface elevation changes during 
the periods of 2001—2006 and 2006-2012 was established. 

3.2.4 Solar radiation modelling 

The SWin is an important factor affecting glacier melting. The spatial distribution of the SWin on 
the surface of a glacier may result in differences in glacier melting. Therefore, we used a solar 
radiation model to explore the spatial variation of the SWin in Urumqi Glacier No. 1 (Kumar et al., 
1997) and analyzed the SWin of the east and west branches. The solar radiation model included the 
effects of terrain shading and cloud cover. It can compute clear-sky direct and diffuse shortwave 
radiation in response to the geographical position, altitude, slope, and aspect of the glacier (Kumar 
et al., 1997). The forcing data comprised latitude, spatial resolution of the DEM, and ground 
reflectance, with a default of 0.2. The potential SWin was computed and then corrected for cloud 
cover through the SWin measured at an automatic weather station (4025 m a.s.l.) deployed on 
Urumqi Glacier No. 1 from 30 April 2018 to 30 April 2019. The rix,y is the ratio of the potential 
SWin to the measured SWin. The spatially distributed SWin (Qcorrx, y) can be expressed as the 
measured SWin multiplied rix, y (the ratio of the potential SWin to the measured SWin). In our study, 
the UAV-DEM was derived from UAV in 2018 (see Section 3.1.3). To improve the efficiency of 
the computer operation, we resampled the UAV-DEM to a resolution of 30 m, and then carried out 
coordinate normalization using the WGS84 datum (EPSG: 4326) and UTM 45N projection and 
applied the data to obtain the spatially distributed SWin. 

3.2.5 Statistical analysis 

In this study, the statistical analysis of the annual average temperature and annual precipitation 
was carried out using the linear fitting method. At the same time, Origin 2018 (OriginLab, 
Northampton, United States) and ArcMap v.10.3 (Environment System Research Institute, ESRI, 
Redlands, USA) were used to draw pictures in this study. 


4 Results 


4.1 Changes in the two branches 


Figure 3a shows the annual mass balance changes in the east and west branches from 1992/1993 
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to 2018/2019. The results indicate that after the separation of Urumqi Glacier No. 1, the mass 


balance of the east branch became more negative than that of the west branch from 1992/1993 to 


2018/2019, except for 1996/1997 and 2003/2004; both the east and west branches experienced an 


overall decreasing trend in their mass balance with a large inter-annual variability, but the west 


branch changed less than the east branch. The period 1997 


—2008 seemed to be a stage of 


continuous decrease, with an abnormally low value occurring between 2009 and 2010; 


afterwards, a four 


year increase stage occurred and then a relatively steady stage was followed 
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Fig. 3 Changes in the annual mass balance (a) and maximum ice velocity (b) for the east and west branches of 


Urumqi Glacier No. 1 from 1992/1993 to 2018/2019, as well as the glacier terminus change (c) from 1979/1980 


to 2019/2020 
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after 2015, but the highest increase occurred in 2019 again. We checked the meteorological data 
in the years with abnormally high and low values of mass balance and found that summer 
(June-August) temperature was the most important factor, followed by annual average 
temperature and annual precipitation, which were consistent with the discussion of a previous 
study on the relative importance of meteorological parameters (Wang et al., 2016). For instance, 
summer temperature was usually higher than 4.2°C during the study period (except for 2009, 
2014, and 2019, with the values of 3.0°C, 3.2°C, and 3.3°C, respectively), corresponding to 
extremely high mass balance values in these years. 

The maximum ice velocity series from 1992/1993 to 2018/2019 can be seen in Figure 3b, 
which shows how the ice velocity changed after the separation of Urumqi Glacier No. 1, as well 
as the differences between the west and east branches. The patterns of the average and minimum 
ice velocities were similar with the maximum ice velocity. The ice velocities of both the east and 
west branches were in a continuously decreasing state, with the decreasing rate of the west branch 
being larger than that of the east branch overall until 2010/2011, after which it became almost the 
same for the two branches. The higher decreasing rate of the ice velocity in the west branch was 
attributed to the higher slope (Wang et al., 2016). 

Figure 3c shows the changes in the location of the glacier terminus position since 1980. 
Generally, the glacier terminus retreat was larger after the separation of Urumqi Glacier No. 1 
than before. The glacier terminus generally retreated at a rate of 4.0 m/a from 1980 to 1993 
(before the separation), and the retreat rate of the east and west branches was 5.0 and 5.8 m/a, 
respectively, from 1994 to 2020 (after the separation). The glacier terminus change process can be 
roughly divided into three stages: the period before 1993, the period of 1994-2011, and the period 
after 2011. In the first stage (the period before 1993), the glacier terminus retreat rate remained 
steady and low, with a value of approximately 3.60 m/a from 1983 to 1988, while it fluctuated 
significantly in other years. In the second stage (the period of 1994—2011), the east branch had a 
lower glacier terminus retreat rate between 3.1 and 3.7 m/a, whereas the west branch had a much 
higher retreat rate with a stepped fluctuation between 4.7 and 7.0 m/a. In the third stage (the 
period after 2011), the pattern of the glacier terminus retreat rate in the east and west branches 
was reversed compared to the previous stage. The glacier terminus retreat rate in the east branch 
speeded up suddenly in 2011/2012 and remained at nearly 8.0 m/a in other years; in contrast, the 
glacier terminus retreat rate in the west branch slowed down to less than 4.0 m/a. In recent years, 
the glacier terminus retreat has become increasingly complex. 

The ice thickness change trends in different periods are shown in Figure 4. The thinning rate of 
the entire glacier was calculated using the area-weighted average of the two branches. 
Specifically, the thinning rate of the east branch was larger than that of the west branch on 
average, but both of them remained almost the same in the early years following the separation of 
the glacier. Since 2006, the thinning rate of the east branch was remained clearly larger than that 
of the west branch, but the difference between the east and west branches seemed to have 
decreased in recent years. 

Photographs of Urumqi Glacier No. 1 in different years (Fig. 5) can visually show the glacier 
shrinkage process of the glacier. Table 2 presents the glacier area change values of the east and 
west branches from 1994 to 2018. These values demonstrate that the area of the east branch 
reduced from 1.115 to 0.983 km? at a rate of 5.740x10° km?/a, while the west branch reduced 
from 0.627 to 0.559 km? at a rate of 2.960x10 km?/a. The relative change rates were —0.006 
km?/a for the east branch and —0.003 km?/a for the west branch, indicating a large difference 
between the two branches. However, the shrinkage of the west branch was larger than that of the 
east branch before 2009 but decreased thereafter. The main reason is that the narrow shape and 
southeast- and south-dominated aspects of the west branch tongue in the earlier period were 
beneficial for glacier terminus melting and retreat. The glacier area shrinkage rate was much 
larger after the separation of the glacier (0.009 km?/a) than that in the previous period from 1962 
to 1992 (0.002 km7/a). 
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Fig. 5 Photographs of Urumqi Glacier No. 1 during different periods showing the glacier shrinkage over the 
past 60 a. (a), 1962; (b), 1993; (c), 2021. 


Table 2 Area changes of Urumqi Glacier No. 1 from 1994 to 2018 


. East branch West branch Urumqi Glacier No. 1 
SEND Area (km?) Change rate (km?/a) Area (km?) Change rate (km7/a) Area (km?) Change rate (km/a) 
1994 1.115 0.627 1.742 
2000 1.111 —0.001 0.622 —0.001 1.733 —0.001 
2001 1.101 —0.010 0.607 -0.015 1.708 —0.025 
2006 1.086 —0.003 0.591 —0.003 1.677 —0.006 
2009 1.069 —0.006 0.578 —0.004 1.647 —0.010 
2012 1.020 —0.016 0.570 —0.003 1.590 —0.019 
2015 0.994 —0.009 0.561 —0.003 1.555 —0.012 
2016 0.991 —0.003 0.559 —0.002 1.550 —0.005 
2017 0.984 —0.007 0.558 —0.001 1.542 —0.008 
2018 0.970 —0.014 0.551 —0.007 1.521 —0.021 
1994-2018 —0.006 —0.003 —0.009 


4.2 Topographic differences of the two branches 


Combined with glacier area, terminus (or length), hypsometry, and accumulation area ratio, the 
topographic characteristics of the east and west branches could be elucidated by their maps of the 
altitude, slope, aspect, and subglacial topography, as shown in Figures 6 and 7. The altitude of 
Urumqi Glacier No. 1 ranged from 3740 to 4478 m a.s.l., where the average altitudes of the east 
and west branches were 3994 and 4127 m a.s.l. in 2018, respectively. In comparison, the average 
altitudes were 3985 and 4123 m a.s.l. for the east and west branches in 1994, respectively. Both 
the average and terminus altitudes of the east branch were lower than those of the west branch. 
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The average slope angles of the east and west branches were approximately 22° and 24° in 2018, 
and about 21° and 22°, respectively, in 1994. In contrast, both the slope distribution and average 
slope value show that the slope angle of the east branch was lower than that of the west branch. 
The aspect maps show that the west branch faced the east, whereas the east branch mainly faced 
the north. In addition, the subglacial topography, as shown in Figure 7, indicates that the basal 
topography of the east branch was more undulating than that of the west branch and that the ice 
thickness of the east branch was larger than that of the west branch. 

Our results further demonstrate that the shrinkage of the west branch was larger than that of the 
east branch before 2009, and the narrow shape and southeast- and south-dominated aspects of the 
west branch tongue during the same period were beneficial for glacier terminus melting and 
retreat; therefore, the slope and aspect of the glacier strongly influenced the glacier area 
reduction. In addition, the shape of the glacier directly affected the ice dynamics. Urumqi Glacier 
No. 1 exhibited wide cross-sections in the firn basins and narrow cross-sections in their ablation 
zones (Fig. 6). The ice velocity, an ice dynamic element, was also related to the ice thickness, and 
the distribution pattern of the ice thickness was similar with that of the ice velocity (Fig. 3) (Wang 
et al., 2018). 
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Fig. 6 Altitude (a and d), slope (b and e), and aspect (c and f) in the east and west branches of Urumqi Glacier 
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Fig. 7 Surface elevation and subglacial topography for the main lines of the east branch (a) and west branch (b) 
of Urumqi Glacier No. 1 
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4.3 Impact of the incoming shortwave radiation (SWi,) on the glacier changes 


During the period from 30 April 2018 to 30 April 2019, the spatial variability of the SWin across 
the glacier surface became much larger as the altitude increases, based on the VAV-DEM shown 
in Figure 8. The mean change of the SWin calculated by the UAV-DEM for all combined 
topographic components was the greatest at lower altitudes (<3886 m a.s.l.) and then began to 
decrease with the increase of altitudes; however, at higher altitudes (>4300 m a.s.1.), it exhibited a 
slightly increasing trend. The mean changes of the SWin in the west and east branches were 143 
and 150 W/m?, respectively. The SWin in the east branch was slightly larger than that in the west 
branch, which may be due to the low altitude of the east branch (Fig. 9). Additionally, the 
estimated result of the UAV-DEM indicates that the northern aspect could receive more SWin, as 
shown in Figures 6f and 8. Such variations would surely have a significant effect on the heat 
budget of different sites, thus influencing the glacier surface melting. 
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Fig. 8 Spatial distribution of the mean changes in the incoming shortwave radiation (SWin) in the east and west 
branches of Urumqi Glacier No. 1 from 30 April 2018 to 30 April 2019 


5 Discussion 


5.1 Uncertainties of individual glacier parameters 


The glacier mass balance within a budget year can be obtained by using the stakes or snow pits. In 
this study, we aimed to analyze the difference in the mass balance of the two branches before and 
after 1993 (the separation of Urumqi Glacier No. 1) to potentially reflect the impact of 
topography on the glacier changes. Uncertainties of the in situ mass balance can usually be 
categorized into three terms. First, according to Fountain and Vecchia (1999) and Kuhn (1999), a 
value of +0.10 m w.e./a is usually adopted as field measurement. Second, Hock and Jensen (1999) 
demonstrated that the error introduced by the Kriging interpolation method is approximately 
+0.10 m w.e./a for a single point mass balance. Third, the glacier retreat could result in an 
overestimation of mass loss of +0.02 m w.e./a for the Urumqi Glacier No. 1 (Xu et al., 2019). 
Additionally, internal accumulation and ablation cannot be quantified. We also recognized that 
internal accumulation and ablation could lead to small but significant systematic cumulative 
errors. By considering the factors mentioned above, the annual uncertainty of the in situ mass 


ChinaX iva (E#AF! 


LI Hongliang et al.: Effect of topography on the changes of Urumqi Glacier... 


250 
= —— East branch 
E£ 2004 
Z 
S | | | | 
Z 1504 | MN E VAN 
o 
=] 
=} 
o 1004 
èn 
5 
a 
[>] 
5 504 
O 
= 
0 T T T T T T T 
3800 3900 4000 4100 4200 4300 4400 
250 
— West branch 
= 2004 
a 
z 150 4 
5 
a 
o 1004 
ob 
S 
a 
[5] 
g 50 4 
= 
0 T T T T T T T 
3800 3900 4000 4100 4200 4300 4400 


Altitude (m) 


Fig. 9 Mean changes in the SWin of the east branch (a) and west branch (b) of Urumqi Glacier No. 1 along with 
the altitude from 30 April 2018 to 30 April 2019 


balance was therefore +0.14 m w.e./a without taking into account the systematic uncertainties. 
The cumulative uncertainty of the mass balance was +0.73 m w.e. during the period of 
1993-2019. 

To determine the glacier surface elevation changes over the respective time periods, we 
compared the DEMs of the two branches at different periods, as mentioned in Section 3.1.3. 
According to Wang et al. (2014), the relative error of the DEMs during the period of 1994—2009 
was less than +0.30 m based on the GPS points in the non-glaciated area. For the DEMs between 
2012 and 2015, ten independent GPS points in the non-glaciated area were selected to perform the 
uncertainty of the DEMs with +0.33 m, and after the co-registration, the mean surface elevation 
difference in the non-glaciated area was +0.41 m (Li et al., 2022). For the DEMs between 2015 
and 2017, the standard error was estimated by calculating the surface elevation changes from the 
terrestrial laser scanning system over the non-glaciated area with 0.23 m in 2015-2016 and 0.16 m 
in 2016-2017 (Xu et al., 2019). The surface elevation difference in the non-glaciated area during 
2017-2018 was 0.68 m. It was concluded that the surface elevation difference in the 
non-glaciated area had stabilized, making the processed DEM difference suitable for estimating 
the glacier surface elevation changes. 

The glacier terminus position data were input into the same coordinate system, and the glacier 
terminus change was obtained by comparing data from various periods. The tape measurement 
error was very small and mainly determined by manual operation. The area change of the glacier 
can be derived from a comparison of the topographic maps for different periods. After the 
integrated evaluation, the whole uncertainty in the glacier area change during the period of 
1994-2018 was approximately 0.01 km? according to Williams et al. (1997), Hall et al. (2003), 


JOURNAL OF ARID LAND 


and Silverio and Jaquet (2005). 

The accuracy of the ice thickness estimation depended on the accuracy of the measurement 
system and the properties of the ice and bedrock. The relative error of the measurement system 
was approximately 1.18%. Combined with the survey points and measurement path differences 
for the different surveys, the ice thickness distribution for each survey caused errors by using 
different interpolation methods. The properties of the ice and bedrock (e.g., crevasses and bedrock 
roughness) had a significant influence on ice thickness estimation. However, it is difficult to 
determine this. The error caused by neglecting the firn layers was estimated at approximately 
5.00%, according to Haeberli et al. (1982). 

The relationship between the ice thickness and surface elevation changes was compared in this 
section. The comparison results show that the correlation coefficient was 0.75 and the root mean 
square error was 3.40 m. The difference was within 10.00%, indicating a good agreement between 
the ice thickness and surface elevation changes. Therefore, the ice thickness change trend can be 
demonstrated in more detailed time intervals over the long-term period. 

The estimated error of the ice velocity at each point was <10.00% of the ice velocity 
components in the X and Y axis directions (Wang et al., 2018). The uncertainty of the ice velocity 
and its changes also depended on the interpolation method and glacier surface characteristics, 
such as the ice thickness and slope. 


5.2 Effect of topography on the individual glacier parameters 


5.2.1 Glacier mass balance and ice thickness 

If adjacent glaciers or branches of a glacier have the same topography, their mass balance should 
be close because they are in the same regional climate conditions. However, differences in the 
topography (see Fig. 6) can result in the differences in the mass balance between them. Affected 
by the altitude, slope and aspect of the glacier, the energy received by the glacier is different, 
which eventually leads to the difference in the glacier mass balance. 

First, because the average altitude of the east branch is lower than that of the west branch, the 
increase of the SWin will inevitably lead to a greater amount of glacier melting on the surface of 
the east branch than on the surface of the west branch (Fig. 9). Second, we calculated the 
accumulation area ratio of the two branches annually and found that the accumulation area ratio 
of the east branch was generally lower than that of the west branch, as shown in Figure 10a; 
therefore, the larger ablation area ratio of the east branch resulted in a higher glacier mass loss 
rate. Third, both east and south aspects influenced the amount of solar radiation received. Figure 6 
shows that the lower section of the west branch faced the east and south more than the east 
branch, suggesting more intensive glacier melting in the west branch than in the east branch. 
However, a large part of the upper section in the east branch facing the northwest was beneficial 
for receiving strong solar radiation in the afternoon, whereas the upper section of the west branch 
faced mainly the northeast and north, which means that it received less solar radiation (Fig. 9). 
According to the study of Yue et al. (2017), on average, the albedo did not have an apparent 
difference between the east and west branches, although it was slightly less for the west branch. 
Therefore, it can be concluded that the lower altitude and accumulation area ratio are the main 
factors leading to a more negative glacier mass balance in the east branch. The influence of 
topography on the solar radiation is complex due to the complicated surface features and aspect 
changes in the lower section of the branches over time. For example, the area of the west branch 
tongue facing the south shrank as the glacier terminus retreat. The albedo was closely related to 
both the surface impurity concentration and topographic shield. 

The topographic difference between the east and west branches led to a similar trend in the ice 
thickness and glacier mass balance changes, that is, the thinning rate of the east branch was larger 
than that of the west branch on average. This is reasonable because the ice thickness change is 
determined mostly by the glacier mass balance (Huang et al., 2018), so a higher glacier mass loss 
will cause a greater thinning in the east branch than in the west branch. 


ChinaXivA ERAT! 


LI Hongliang et al.: Effect of topography on the changes of Urumqi Glacier... 


—— East branch 


—e— West branch 


(a) 


100 


T 
© 
a 


- — &, 
seceooeooe eco oc G7 
oreonrAnana 


(%) onel eaw uoyemumooy 


—— East branch 
—e— West branch 


(b) 


S 
T 


- 

n 

= 
l 


-2.0 | 
-3.0 | 


N 
T 


(B/W) 9781 ANAI SNUTULID) IAMJ O 


East branch 


—e— West branch 


— 


0 
3.5 


S 
oa) 


e n Se n 
N = 


n 
N = S 


(eu) SNUTULID) 131923 OY} JO AWIOOJOA 9d] 


LAGE LOO’ Z0cc0d-AIXeUIYS 


ar 


e 


Fig. 10 Accumulation area ratio (a), glacier terminus retreat rate (b), and ice velocity of the glacier terminus (c) 


in the east and west branches of Urumqi Glacier No. 1 from 1993/1994 to 2017/2018 


5.2.2 Ice velocity 


Because the accumulated glacier mass in the upper section of the two branches is transported 
down to the ablation area through ice movement, the ice velocity is one of the significant factors 


controlling the surface and dimensional changes of Urumqi Glacier No. 1. The slope influences 


the ice velocity and mass balance, and may also impact the snow accumulation rate through its 
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effect on the avalanche transport of snow on the glacier surface. Figure 3b shows that the ice 
velocity of the west branch was always higher than that of the east branch, but the difference 
decreased, particularly after 2010. Previous studies, such as Li et al. (2011), concluded that 
Urumqi Glacier No. 1 was predominantly cold-based, with basal sliding occurring only close to 
its snout. It is well known that the ice velocity caused by creep deformation mainly depends on 
the ice thickness and slope (Paterson, 1994). After analyzing the effects of the ice thickness and 
slope using the basic dynamic principle, Wang et al. (2018) found that the ice velocity is more 
sensitive to the slope compared to the ice thickness, and a larger slope will lead to a higher ice 
velocity in the west branch than in the east branch. Figure 6 shows that the slope was higher in the 
west branch than in the east branch. The calculated average slope indicates that the difference in 
the slope between the west and east branches declined to some extent, as mentioned above. This 
coincides with the change trend of the maximum ice velocity, as shown in Figure 3b; that is, the 
ice velocity was higher in the west branch than in the east branch, gradually becoming closer 
under current slope changes. Because the average ice thickness in the early 2000s was larger in 
the east branch (approximately 59.0 m) than in the west branch (approximately 45.0 m), the 
contribution of the ice thickness to the ice velocity was larger in the east branch than in the west 
branch. However, the measured ice velocity in the east branch was lower than that in the west 
branch, indicating that the effect of the ice thickness was less than that of the slope. 

From Figure 7 we can conclude that compared with the west branch, the more undulant basal 
topography of the east branch has a retardative effect on the ice movement. Therefore, the 
difference in the ice velocity between the west and east branches is a result of the comprehensive 
effects of the slope, ice thickness, and basal topography. 


5.2.3 Glacier terminus 


Changes in the terminus of a glacier are mainly controlled by the glacier mass balance and ice 
velocity near the terminus. The results of this study mentioned above indicate that the glacier 
mass balance of the east branch was more negative, and the ice velocity of the east branch was 
lower than that of the west branch, which makes the mass supply in the terminus of the east 
branch slow. The combination of the above two factors led to a larger glacier terminus retreat for 
the east branch than for the west branch. However, before 2012, the measured glacier terminus 
retreat was lower for the east branch than for the west branch (Wang et al., 2016; Li et al., 2022). 
Comparing the ice thickness and shape of the east branch with the west branch, the west branch 
tongue was much narrower with a pointed terminus, making the absolute ice thickness thinner 
than that of the east branch. Another factor is that the southeast aspect of the west branch tongue 
in the earlier period benefitted more from the solar radiation for glacier melting; therefore, the 
glacier melting rate of the west branch was stronger (Wang et al., 2016). This feature must have 
led to a larger terminus retreat rate in the west branch compared with the east branch if the 
differences in the glacier melting rate and ice velocity between the east and west branches were 
not too large. From the observation data in Figure 10b and c, the glacier melting rate did not show 
any system differences at the stakes closest to the terminus of the east or west branch, while the 
ice velocity at the stake of the west branch was much higher than that of the east branch before 
2007, but the difference gradually decreased thereafter. Thus, the glacier terminus retreat was 
influenced by the combined effects of the tongue shape, melting rate, and ice velocity. 

Since 2011/2012, the difference of the glacier terminus retreat between the west and east 
branches seemed to have become more complex. The west branch retreated more than the east 
branch in some years but less in other years, with frequent alternation after 2015. In view of the 
small difference in the ice velocity between the west and east branches after 2010, the glacier 
melting rate and topography were the main causes of the difference in the glacier terminus retreat. 
In particular, the difference of glacier melting was more important because the topography 
changed little from year to year. 


5.3 Changes of the east and west branches response to climate change 


It is widely accepted that glaciers in the High Mountain Asia will continue to shrink under the 
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background of climate warming or cessation of warming due to the lag of glacier responses to 
climate change (Li et al., 2011; Yao et al., 2012; Farinotti et al., 2015; Zemp et al., 2015, 2019; 
IPCC, 2019), as is the case for Urumqi Glacier No. 1 (Li et al., 2011; Xu et al., 2011; Wang et al., 
2016). However, it still needs to explore whether this glacier retreat will accelerate, although 
some studies have indicated that the glacier retreat in the Tianshan Mountains, including or 
specifically Urumqi Glacier No. 1, has accelerated in recent decades and will continue to do so in 
the near future (Pieczonka et al., 2013, 2015; Brun et al., 2017). The observational data suggest 
that this glacier shrinkage became more rapid after the 1990s than before, but the observed glacier 
mass balance did not appear to have become more negative since 2012, based on a comparison 
between the mean values of 2012-2019 and 1997-2011 (Fig. 3a). The average glacier mass 
balance was —0.8 m w.e. during the period of 1997-2011 and —0.7 m w.e. during the period of 
2012-2019 in the east branch, whereas it was —0.6 and —0.4 m w.e. during the periods of 
1997-2011 and 2012-2019 in the west branch, respectively. If the current climate conditions are 
maintained in the future, the negative glacier mass balance and glacier terminus retreat will 
undoubtedly continue, but the glacier terminus retreat will slow gradually when the terminus 
positions of the east and west branches retreat at higher altitudes. 

Owing to the topographic differences, changes in the glacier mass balance, ice velocity, and 
glacier terminus retreat will inevitably differ between the east and west branches, despite the same 
climate conditions. As mentioned above, the glacier mass balance was more negative for the east 
branch than for the west branch because of the lower average altitude and accumulation area ratio 
in the east branch. The ice velocity of both the east and west branches will continue to decrease in 
the future, but the velocity of the east branch may catch up and even exceed that of the west 
branch, because the east branch is thicker on average and the slope difference between the west 
and east branches diminishes gradually. The difference of the glacier terminus retreat is 
complicated because it is dependent not only on the glacier melting rate and ice velocity at the 
glacier terminus but also on the shape of the glacier's snout. In addition, the shape and size of the 
glacier are important factors in controlling the glacier dynamics. It was previously found that 
small glaciers lost more area than large glaciers, because small glaciers such as Urumqi Glacier 
No. | are generally more sensitive to climate change (Oerlemans, 1989). The response time of the 
glacier to climate change can be roughly estimated from the simplest formula introduced by 
Paterson (1994), which is based on an ideal assumption. According to Wang et al. (2016), the 
actual response time is likely to be shorter than the estimated value, perhaps around 10 a. This 
suggests that Urumqi Glacier No. 1 will continue to shrink in the coming decade, even without 
further climate warming, due to the time it takes to respond to the previous climate warming. 
Nevertheless, it can be imagined that when the east and west branches shrink, the glacier terminus 
retreat of the east branch will be faster than that of the west branch, as the lower altitude and 
western-facing direction of the east branch would lead to more ice loss in comparison to the west 
branch. 


6 Conclusions 


Under the background of climate warming, Urumqi Glacier No. 1 has experienced continuous 
mass loss and dimensional shrinkage in recent decades; however, after the separation of the 
glacier into two branches, no obvious accelerated trend was identified in recent years. The 
observation data show that there were differences between the east and west branches in terms of 
glacier mass balance and dimensional changes, such as terminus retreat, thickness decrease, and 
ice velocity change. The difference of the glacier area shrinkage between the two branches was 
not large, although the east branch shrunk slightly more. From the comparison of the topographic 
features of the two branches, the greater negative glacier mass balance in the east branch was 
mainly attributed to the lower average altitude and accumulation area ratio. As expected, the 
larger negative glacier mass balance has led to a higher average thinning rate of the east branch, 
although the thinning rate between the east and west branches differed in different years. The 
narrow and pointed snout of the west branch and its southeast aspect resulted in a larger glacier 
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terminus retreat for the west branch. The ice velocity decreased in both the east and west 
branches, but the decrease rate was lower in the east branch owing to its increase in the slope. 

A topographic solar radiation model was used to quantify the spatial variability of the SWin 
across the glacier surface. It is indicated that the SWin became much larger as altitude increased 
based on the UAV-DEM, and the north aspect could receive more SWin, resulting in glacier 
melting. Importantly, the SWin was the greatest at lower altitudes (<3886 m a.s.l.) and then began 
to decrease with the increase of altitudes; however, at higher altitudes (>4300 m a.s.l.), it 
exhibited a slightly increasing trend. The SWin received by the east branch was slightly larger 
than that received by the west branch, which was probably due to the lower altitudes of the east 
branch. 

In the future, both branches of Urumqi Glacier No. 1 will continue to shrink due to the climate 
warming; however, differences between the two branches will still exist owing to their 
topographic differences. The east branch retained a more negative glacier mass balance and a 
larger thinning rate than the west branch. The ice velocity of both the east and west branches will 
decrease continuously, but the velocity of the east branch may catch up to and even exceed that of 
the west branch after some years. With the continuous decrease in the glacier area, the glacier 
terminus retreat of the east branch may exceed that of the west branch, because the east branch 
has a lower average altitude and the western-facing direction of the east branch would lead to 
more ice loss in comparison to the west branch. 
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